pub fn new_birthday_probability(n: u32) -> f64 {
    let mut p = 1f64;
    for i in 0..n {
        p *= (365f64 - i as f64) / 365f64;
    }

    1f64 - p
}
